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Abstract: 

The 16-22 amino acid fragment of the /3-amyloid peptide associated with the Alz- 
heimer's disease, A/3, is capable of forming amyloid fibrils. Here we study the ag- 
gregation mechanism of A/?i6_22 peptides by unbiased thermodynamic simulations at 
the atomic level for systems of one, three and six A/?i6-22 peptides. We find that 
the isolated A/5i6_22 peptide is mainly a random coil in the sense that both the a- 
helix and /3-strand contents are low, whereas the three- and six-chain systems form 
aggregated structures with a high /5-sheet content. Furthermore, in agreement with 
experiments on A/5i6_22 fibrils, we find that large parallel /5-sheets are unlikely to 
form. For the six-chain system, the aggregated structures can have many different 
shapes, but certain particularly stable shapes can be identified. 



*E-mail: favrin, anders, sandipan@thep.lu.se 
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1 Introduction 



The fibrillar aggregates that characterize amyloid diseases, such as the Alzheimer's 
disease, are formed by specific peptides or proteins. However, it is known that several 
non-disease- related proteins are capable of forming similar amyloid structures [1,2], 
and that the aggregation of such proteins can be cytotoxic [3]. This suggests, first, that 
polypeptide chains have a general tendency to form amyloid structures and, second, 
that natural proteins should have evolved mechanisms to avoid this tendency. Such 
mechanisms have indeed been proposed [4-6] . The propensity of a given polypeptide 
chain to form amyloid fibrils depends, nevertheless, on its amino acid sequence [7-11], 
and short sequence stretches promoting amyloid formation have been identified [12, 
13]. 

While the structure of amyloid fibrils is not known in atomic detail, there is ample 
evidence from X-ray fiber diffraction studies that the core of the typical amyloid fibril 
is composed of /9-sheets whose strands run perpendicular to the fibril axis [14]. More 
detailed information is available, for example, for fibrils made from difi^erent fragments 
of the Alzheimer's A/3 peptide. In particular, there is evidence from solid-state NMR 
studies for a parallel organization of the /3-strands in A/?io-35 [15] and A/3i_4o [16] 
fibrils, and for an antiparallel organization in A(3s4-4:2 [17], A(3u-25 [18] and A/3i6-22 
fibrils [19,20]. Most of these fragments contain the hydrophobic A/3i6-20 segment 
(KLVFF), which is known to be important in the A/S-A^ interaction [21]. 

Small peptides hke A/3i6-22 are well suited as model systems for probing the mecha- 
nisms of aggregation and fibril formation, and are being studied not only in vitro but 
also in silico. Computer simulations of simplified [22-26] and atomic [27-31] models 
have provided useful insights into the aggregation behavior of some peptide systems. 
To properly explore the free-energy landscape of aggregation at the atomic level is, 
nevertheless, a computational challenge. 

Here we investigate the formation and properties of A/9i6-22 oligomers by unbiased 
Monte Carlo (MC) simulations of systems with up to six chains, using a sequence- 
based atomic model with an effective potential based on hydrogen bonds and hy- 
drophobic attraction (no explicit water molecules). The same model has previously 
been used to study the folding of individual peptides [32-34] . It was shown that this 
model is able to fold several different peptides, both a-helical and /?-shcct peptides, 
for one and the same choice of parameters. The calculated melting behaviors were, 
moreover, in good agreement with experimental data for all these peptides. 



2 



2 Model and Methods 



The main object of study in this paper is the peptide A/?i6_22, given by acetyl-Lys- 
Leu-Val-Phe-Phe-Ala-Glu-NH2. We consider systems of one, three and six A/?i6-22 
peptides. The multichain systems are contained in periodic boxes. All the interactions 
are short range, which makes the implementation of the periodic boundary conditions 
straightforward. The box sizes are (35A)^ and (44A)^ for three and six chains, re- 
spectively, corresponding to a constant peptide concentration. For computational 
efficiency, the peptide concentration is taken to be high. 

Our model [32-34] contains all atoms of the peptide chains, including hydrogen atoms. 
The model assumes fixed bond lengths, bond angles and peptide torsion angles (180°), 
so that each amino acid only has the Ramachandran torsion angles 0, ip and a num- 
ber of side-chain torsion angles as its degrees of freedom. Numerical values of the 
geometrical parameters held constant can be found elsewhere [32]. 



The interaction potential 

E = + E\oc + -Ehb + -£^hp (1) 

is composed of four terms, which we describe next. Energy parameters are given 
on a scale such that a temperature of T = 300 K corresponds to kT ^ 0.447 {k is 
Boltzmann's constant) [34]. 



The first term in Eq. ^ E^^, represents excluded- volume effects and has the form 

£^ev = ft^ev [ ''^^ ^ — , (2) 

where the summation is over pairs of atoms (i, j), ^ev = 0.10, and cTj = 1.77, 1.75, 
1.55, 1.42 and 1.00 A for S, C, N, O and H atoms, respectively. The parameter \ij has 
the value 0.75 for all pairs except those connected by three covalent bonds, for which 
\ij = 1. When the two atoms belong to different chains, we always use Xij = 0.75. 
To speed up the calculations, Eq. |21is evaluated using a cutoff of r^j = 4.3Ajj A, and 
pairs with fixed separation are omitted. 



The second energy term, E\oc, is a local intrachain potential. It has the form 
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where the inner sum represents the interactions between the partial charges of the 
backbone NH and CO groups in one amino acid, /. This potential is not used for Gly 
and Pro amino acids which have very different 0, ip distributions, but is the same for 
all other amino acids. The inner sum has four terms (NC, NO, HC and HO) which 
depend only on the and if) angles for amino acid /. The partial charges are taken 
as qi = ±0.20 for H and N and = ±0.42 for C and O [35], and we put kiqc = 100, 
corresponding to a dielectric constant of ~ 2.5. 

The third term of the energy function is the hydrogen-bond energy E^y,, which has 
the form 

^hb = u{rij)v{aij,(3ij) + eg] ^ u{rij)v{aij, , (4) 

bb— bb sc— bb 

where the two functions u{r) and v(a,P) are given by 

, f (cos a COS ifa,/3>90° 

= \0 otherwise 

We consider only hydrogen bonds between NH and CO groups, and Vij denotes the 
HO distance, aij the NHO angle and (3ij the HOC angle. The parameters e^^, e^^^ 
and (Jhb are taken as 3.1, 2.0 and 2.0 A, respectively. The function u{r) is calculated 
using a cutoff of = 4.5 A. The first sum in Eq. |3] contains backbone-backbone 
interactions, while the second sum contains interactions between charged side chains 
(Asp, Glu, Lys and Arg) and the backbones. For intrachain hydrogen bonds we make 
two restrictions. First, we disallow backbone NH (CO) groups to make hydrogen 
bonds with the two nearest backbone CO (NH) groups on each side of them. Second, 
we forbid hydrogen bonds between the side chain of one amino acid with the nearest 
donor or acceptor on either side of its C^. For interchain hydrogen bonds, we make 
no such restrictions. As a simple form of context dependence, we assign a reduced 
strength to hydrogen bonds involving chain ends, which tend to be exposed to water. 
Following the experimental studies of the A/3i6_22 peptide [19,20], we have used acetyl 
and amide capping groups at the ends. A hydrogen bond involving one or two such 
groups is reduced in strength by factors of 2 and 4, respectively. 

The fourth energy term, i?hp, represents an effective hydrophobic attraction between 
nonpolar side chains. It has the pair-wise additive form 

E^^ = -Y,MuCu, (7) 

KJ 
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I II III 



I Ala 0.0 0.1 0.1 

II He, Leu, Met, Pro, Val 0.9 2.8 

III Phe, Trp, Tyr 3.2 

Table 1: The hydrophobicity matrix M/j. Hydrophobic amino acids are divided into 
three categories. The matrix M/j represents the size of hydrophobicity interaction 
when an amino acid of type / is in contact with an amino acid of type J. 



where C/j is a measure of the degree of contact between side chains / and J, and Mjj 
sets the energy that a pair in full contact gets. The matrix Mjj is defined in Table ^ 
To calculate C/j we use a predetermined set of atoms, A/, for each side chain /. We 
define Cu as 

where the function f{x) is given by f{x) = 1 ii x < A, f{x) = if x > B, and 
fix) = (B - x)/{B - A) if A < X < B [A = (3.5 A)^ and B = (4.5 A)^]. Roughly 
speaking, Cjj is the fraction of atoms in Aj or Aj that are in contact with some 
atom from the other side chain. For Pro, the set Aj consists of the C/3, and 
atoms. The definition of Aj for the other hydrophobic side chains has been given 
elsewhere [32]. For pairs that are nearest or next-nearest neighbors along the same 
chain, we use a reduced strength for the hydrophobic attraction; Mjj is reduced by a 
factor of 2 for next-nearest neighbors, and taken to be for nearest neighbors. 

To study the thermodynamic behavior of this model, we use simulated tempering [36- 
38] in which the temperature is a dynamical variable. For a review of simulated 
tempering and other generalized-ensemble techniques for protein folding, see Ref. [39]. 
We study the one- and three-chain systems at eight different temperatures, ranging 
from 275 K to 369 K, and the six-chain system at seven temperatures, ranging from 
287 K to 369 K. 



Our simulations are carried out using two different elementary moves for the backbone 
degrees of freedom: first, the highly non-local pivot move in which a single backbone 
torsion angle is turned; and second, a semi-local method [40] that works with up 
to eight adjacent backbone degrees of freedom, which are turned in a coordinated 
manner. Side-chain angles are updated one by one. In addition to these updates, we 
also use rigid-body translations and rotations of whole chains. Every update involves 
a Metropolis accept/reject step, thus ensuring detailed balance. All our simulations 
are started from random configurations. All statistical errors quoted are la errors 



5 



Hi 0| H| 0| H 0| 

id ol nI cl n] cJ 
h! o' 




Lys c 



,0 |H |0 |H 

c In ^ I c 



Glu 



O |H 



lO |H 

Jn 



O H 



LysH| O 



O 'H 



Lys 



N| C 

h! o 



H| 0| 
CT NY C~ 

h! o! h! O. 



H| 0| Glu 



.c In ^ I c 



O |H 



C 
H O 



lO iH 

/N 



TC 

Glu Iq 'h 



'o 'h 'o 'h Lys 



Figure 1: Schematic illustrations of the hydrogen bond patterns for in-register, parallel 
/3-strands (left) and in-register, antiparallel /^-strands (right). 

obtained from the variation between independent runs. We performed 9 runs with 
10^ elementary MC steps for A^^c = 1, 11 runs with 10^ MC steps for = 3, and 18 
runs with 2 ■ 10^ MC steps for = 6. Each of the = Q runs required about 12 
CPU days on a 1.6 GHz computer. 

To characterize the behavior of these systems, we first determine the secondary struc- 
ture. For a chain with amino acids, we define the a-helix and /?-strand contents 
as the fractions of the N — 2 inner amino acids with their (0, ip) pair in the a-helix 
and /5-strand regions of the Ramachandran space. We assume that a-helix corre- 
sponds to —90° < (p < —30°, —77° < tp < —17° and that /3-strand corresponds to 
— 150° < < —90°, 90° < ip < 150°. The average a- helix and /5-strand contents, over 
all the chains of the system, are denoted by H and 5, respectively. 

To distinguish between parallel and antiparallel /?-sheet structure, we examine the 
orientation of end-to-end vectors. For a given multichain configuration, we first de- 
termine all pairs of chains such that (i) their interchain hydrogen bond energy is less 
than — 1.5e[jJ,'' (roughly corresponding to 2-3 hydrogen bonds), and (ii) both chains 
have a /9-strand content higher than 0.5. For each such pair of chains, we then cal- 
culate the scalar product of their normalized end-to-end unit vectors. If this scalar 
product is greater than 0.7 (less than —0.7), we say that the two chains are parallel 
(antiparallel). We denote the numbers of parallel and antiparallel pairs of chains by 
n+ and n_, respectively. Fig. illustrates the hydrogen-bond patterns in parallel and 
antiparallel /3-sheets. 
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Figure 2: Monte Carlo evolution in a simulated-tempering run for Nc = Q A/3i6_22 
peptides, (a) The total energy E (full line) and the hydrogen-bond energy i?hb (dashed 
line), both in kcal/mol. (b) The temperature index k. There are seven allowed 
temperatures Tk, satisfying To = 287 K < Ti < . . . < Tq = 369 K. Measurements are 
taken every 10^ MC steps. 



3 Results and Discussion 



Using the model described in the previous section, we study the thermodynamics of 
systems of Nc A/3i6-22 peptides for iVc = 1, 3 and 6. Fig. |21 illustrates the Monte 
Carlo evolution in one of 18 independent simulated-tempering runs for the six-chain 
system. In the course of the run, aggregated low-energy structures form and dissolve 
several times. 
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Figure 3: (a) The a-helix content H against temperature T for A/5i6_22 for A^^c = 1 
(o), Nc = ^ (•) and A^'c = 6 (A). Lines joining data points are only a guide for the 
eye. (b) Same for the /5-strand content S. Note that the scales in (a) and (b) are 
different. 

3.1 Secondary Structure 

Fig. El shows the a-helix and /3-strand contents H and 5, as defined in the previous 
section, against temperature for different A*",,. For Ac = 1, we see that both H and S 
are small at all temperatures studied, although H increases with decreasing tempera- 
ture. So, in our model, the A/3i6-22 monomer is mainly a random coil throughout this 
temperature range. The A"c = 3 and Ac = 6 systems show a qualitatively different 
behavior; S increases sharply with decreasing temperature, to values of S" = 0.6 and 
higher, whereas H is very small. These results clearly demonstrate that unless the 
temperature is too high, the three- and six-chain systems self-assemble into ordered 
structures with a high /?-strand content. 

The temperature at which the aggregation sets in depends strongly on the peptide 
concentration, and exploring that dependence is beyond the scope of the present 
study. We note, however, that the /3-sheet formation sets in at a higher temperature 
for Ac = 6 than for A^c = 3. This fact is also reflected in the behavior of the specific 
heat, as shown in Fig. EJ For Ac = 3 and Ac = 6, the specific heat Cy{T) exhibits a 
pronounced peak. As the system size increases from Ac = 3 to Ac = 6, the peak is 
shifted towards higher temperature. Near the peak, the energy distribution is broad 
(data not shown), showing that both aggregated low-energy and unstructured high- 
energy states occur with a significant frequency at these temperatures. 

Our results for Ac = 1 and Ac = 3 can be compared with results from molecular 
dynamics simulations with explicit water by Klimov and Thirumalai [29]. Using 
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Figure 4: Specific heat against temperature T for Ai'c = 1, 3 and 6 A/3i6-22 peptides, 
as obtained using histogram reweighting techniques [41]. The bands are centered 
around the expected values and show statistical la errors. is defined as = 
{N^N)-^d{E)/dT = {NcNkT^y\{E^) - (E)^), where is the number of chains, 
is the number of amino acids per chain, and (O) denotes a Boltzmann average of 
variable O. 

somewhat different definitions of H and S and a temperature of T = 300 K, these 
authors found that H = 0.11 and S = 0.33 for = 1, and H = 0.26 and S = 0.30 
for Nc = 3. Our A^^c = 1 results (see Fig. ^ are in reasonable agreement with theirs, 
given that we use a stricter definition and /3-strands. However, our A^c = 3 results 
disagree with theirs. They obtained a smaller /^-strand content and a larger a-helix 
content compared to their own A^^ = 1 results; whereas we observe a much larger 
/9-strand content for Ag = 3 compared to Ac = 1. 

For the Ac = 3 system, Klimov and Thirumalai [29] furthermore found evidence for 
an obligatory a-helical intermediate. To see whether or not such an intermediate 
exists in our model, we divided the energy axis into bins and calculated the average 
a-helix and /?-strand contents for each bin, at a fixed temperature near the specific 
heat maximum. Fig. shows the resulting a-helix and /5-strand profiles H{E) and 
S{E). We see that the /3-strand content S{E) increases steadily with decreasing 
energy. The a-helix content H[E), on the other hand, does have its global maximum 
at ~ 130kcal/mol. However, the maximum value of H[E) is very small. Hence, we 
find no sign of an obligatory a-helical intermediate in our model. Most of the amino 
acids in a typical configuration at intermediate energies are either random coils or 
/9-strands. 




9 





40 60 80 100 120 140 

E (kcal/mol) 

Figure 5: The a-helix (o) and /3-strand (•) profiles H{E) and S{E) (see the text) for 
the six-chain A/3i6_22 system at T = 325 K. 

3.2 /^-Strand Organization 

As mentioned in the introduction, there exist experimental results [19,20] suggesting 
that the /3-strands in full A/?i6-22 fibrils have an in-register, antiparallel organization. 
To find out whether our systems show a preference for either parallel or antiparallel 
/3-sheets, we consider the joint probability distribution P{n+,n^), where n+ and n_ 
count the numbers of interacting chain pairs with high /5-strand contents that are 
parallel and antiparallel, respectively (see Sec. |21). 

Table 121 shows this distribution for the A^c = 3 system at T = 275 K. For this system, 
the most probable combination of (n+,n_) is (1,1), corresponding to a mixed (3- 
sheet. At the same time, the distribution shows a clear asymmetry. The frequency of 
occurrence for antiparallel /3-sheets with (n_|_, n_) = (0, 2) is a factor of 7 higher than 
that for parallel /5-sheets with (n_|_,n_) = (2,0). 

The corresponding results for A'c = 6, at T = 287 K, are shown in Table El As 
in the = 3 case, we find that a majority of the configurations contain mixed 
/3-sheet structure, n+ and both being nonzero. The asymmetry of the {n^,n^) 
distribution is even more pronounced for = Q than for = 3. In particular, we see 
that large n_ values are much more probable than large n+ values; the combination 
{n+,n_) = (4,0) is, e.g., very unlikely to occur, whereas (n+,n_) = (0,4) does occur 
with a significant frequency. 

Tables El and El show the {n+,n_) distribution at the lowest temperatures studied. 
With increasing temperature, the average n+ and n_ steadily decrease. At the highest 
temperature studied, 369 K, about 99% of the conformations have n+ = n_ = 0, for 
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i 2 

0.17(2) 0.22(3) 0.14(3) 

1 0.13(2) 0.32(6) 

2 0.020(7) 

Table 2: The probability distribution P(?7,+,n_) (see Sec. |21) for A^c = 3 A/3i6_22 
peptides at T = 275 K. P{n^, n_) values smaller than 10~^ are omitted. The numbers 
in parentheses are the statistical errors in the last digits. 



12 3 4 

0.028 (5) 0.059(11) 0.08 (2) 0.06 (2) 0.030(15) 

1 0.038(6) 0.12(2) 0.16(3) 0.10(3) 0.006(3) 

2 0.026(11) 0.11 (5) 0.14(5) 0.004(2) 

3 0.008(5) 0.013(9) 0.015(12) 

Table 3: Same as Table El for Nc = Q A/5i6-22 peptides at T = 287 K. 



A'^c = 3 as well as A'^c = 6. The full (n+, nJ) distribution for both A'c = 3 and A'c = 6 
at all the different temperatures studied can be found as Supplementary Material. 

Although the statistical uncertainties are somewhat large, the results in Tables |21 and El 
show some clear trends. The most striking one is that large rij^ values are strongly 
suppressed, which means that large parallel /^-sheets are very unlikely to form. The 
probability of having large antiparallel /3-sheets is much higher. Compared to purely 
antiparallel /5-sheet structures, it is possible that mixed /?-sheet structures are more 
difficult to extend to large stable structures. To be able to check whether or not this 
is the case, simulations of larger systems are required. 

Why are antiparallel /3-sheets favored over parallel ones? Klimov and Thirumalai [29] 
concluded that A/?i6-22 peptides make antiparallel /3-sheets because of Coulomb in- 
teractions between charged side chains; the two end side chains of the A/3i6-22 peptide 
carry opposite charges, which indeed should make the antiparallel orientation electro- 
statically favorable. However, our model completely ignores Coulomb interactions be- 
tween side-chain charges and still strongly favors the antiparallel organization. Other 
mechanisms than Coulomb interactions between side-chain charges might therefore 
play a significant role, such as the geometry of backbone-backbone hydrogen bonds 
(see Fig. jT]), steric effects, and the precise distribution of hydrophobicity along the 
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Figure 6: Schematic representations of three different registries for an antiparallel pair 
of A/?i6-22 peptides. 

chains. A recent experimental study [20] highhghts the importance of the hydropho- 
bicity distribution. This study showed that the /5-sheet structure of A/3i6-22 fibrils 
can be changed from antiparallel to parallel by adding an octanyl end group to the 
peptide which increases its amphiphilicity. 

To probe the registry of the /3-sheets, we monitored backbone-backbone hydrogen 
bond patterns (see Fig. [T)). Fig. IHl illustrates three possible antiparallel registries: (a) 
17 + ^ 20 - fc, (b) 17 + A; ^ 21 - fc, and (c) 17 + k ^ 22 - k. The 17 + A; ^ 
21 — k registry is the one found in experiments on A/3i6-22 fibrils [19,20], whereas 
experiments on fibrils made from the slightly larger segment A/3ii_25 found evidence 
for the 17 + k ^ 20 — k registry at pH 7.4 and for the 17 + <-> 22 — A; registry at 
pH 2.4 [18]. In our calculations, the 17 + A; 20 — A; and 17 + k ^ 21 — k registries 
occur with high and comparable frequencies. The 17 + A; ^ 22 — A; registry is, by 
contrast, strongly suppressed, which probably is due to hydrophobic effects, although 
steric clashes between the large Phe side chains could play a role, too. As to the 
17 + A; ^ 20 — A; and 17 + A; ^ 21 — A; registries, it would be very interesting to 
see whether their relative frequencies of occurrence depend on {n^, but that will 
require higher statistics than those provided by the present calculations. 



To test our model, we performed simulations similar to those for the A/?i6_22 peptide 
for some other peptides. Some of these peptides, including the polar one studied by 
Diaz-Avalos et al. [42], had a low overall hydrophobicity. We found that the propensity 
to aggregate is much lower for such peptides than for the A/3i6-22 peptide, and a higher 
peptide concentration was required to promote aggregation. These results clearly show 
that in our model, hydrophobic attraction is a major driving force for aggregation. 

As an example of a peptide with a significant hydrophobicity but an uneven dis- 
tribution of it, we studied the peptide acetyl-Lys-Phe-Phe-Ala-Ala-Ala-Glu-NH2, in 
which the two strongly hydrophobic Phe amino acids are asymmetrically placed. For 



3.3 Other Peptides 
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this peptide, we obtained aggregated (3-sheet structures with a predominantly paral- 
lel /3-strand organization, which in particular confirms that our model is capable of 
generating stable parallel /3-sheets. 



3.4 Examples of Low-Energy Structures 

It is known that relatively small assemblies formed early in the aggregation of full- 
length A/3 [43-45] as well as non-disease-related proteins [3] can be toxic, which makes 
it very interesting to study possible oligomer shapes. In addition, such structures 
represent potential seeds for the fibril formation. 

From our simulations, we find that the six-chain A/5i6_22 system does not exhibit a 
single dominating free-energy minimum, but rather a number of more or less degener- 
ate local minima. Fig. [7| shows two snapshots of such minima. The /3-strand content 
is, as noted earlier, high, and the structures shown in Fig. [3 illustrate this property. 

In the simplest class of typical structures observed in our simulations, five of the chains 
form a relatively flat /5-sheet, whereas the remaining chain is a random coil and held 
in contact with the /?-sheet by hydrophobic attraction. Six-stranded /3-sheets also 
occur in the simulations, but with a low frequency, as can be seen from the P{n+, n_) 
distribution in Table IHl Further, for the six-chain system, we observe the emergence of 
new non-trivial structures with no analogs in the three-chain simulations. The second 
structure in Fig. [7| illustrates this. Here stability is achieved by stacking two different, 
three-stranded, /^-sheets together, which brings hydrophobic side chains from the two 
/^-sheets in close contact. Such "sandwiches" occur with a non-negligible frequency 
in our simulations. To estimate the precise populations of these minima is difficult. 
However, five-stranded /3-sheets did occur more frequently than sandwiches in the 
simulations. By visual inspection, we further estimate that of the order of 10 % of 
the configurations are sandwich-like at the lowest temperature studied, at which the 
snapshots were taken. These low-energy structures also occur at higher temperatures, 
but become very rare above the specific heat maximum (see Fig. 0]). 

In none of our simulations did we find any indication of a free-energy minimum in 
which the /3-strands are joined end-to-end to form the so-called /?-helix [47]. In our 
model, stability is enhanced by increasing the number of hydrogen bonds or by in- 
creasing hydrophobic contacts. For system sizes as small as those we examined, the 
(3-helix is inferior to many competing structures in both these respects, and hence its 
absence is expected. 
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Figure 7: Two typical low-energy structures from our simulations of six A/9i6-22 pep- 
tides: a five-stranded /3-sheet (left), and two three-stranded /9-sheets "sandwiching" 
several of their hydrophobic side-chains between them (right). Drawn with Ras- 
Mol [46]. 



4 Conclusion 



Using a sequence-based atomic model which was originally developed for folding stud- 
ies of single peptides [32-34], we studied the aggregation properties of A/9i6-22 pep- 
tides. In this model, we found that A/?i6-22 peptides have a high propensity to self- 
assemble into aggregated structures with a high /3-strand content, while the isolated 
A/9i6-22 peptide is mainly a random coil. Both parallel and antiparallel arrangements 
of the /5-strands occur in the model, with a definite preference for the antiparallel 
arrangement. 

It is important to note that we find this preference for the antiparallel /3-strand orien- 
tation despite ignoring the Coulomb interactions between the two charged side chains 
at the ends of the peptide. It has been suggested [29] that such Coulomb interactions 
are the main determinant for the antiparallel orientation. While these Coulomb in- 
teractions might enhance the tendency for A/5i6-22 peptides to form /3-sheets with an 
antiparallel organization, our results strongly suggest that other factors play a signif- 
icant role, too. It is worth noting that the orientation is not necessarily determined 
solely by sequence-specific side-chain interactions, as antiparallel /5-sheets are widely 
held to be intrinsically more stable than parallel ones. For the A/5i6_22 peptide, which 
in particular lacks a clear amphiphility, there is no obvious mechanism to overcome 
this tendency. 

In our simulations, we did not observe an absolute free-energy minimum, but rather 
several nearly degenerate minima corresponding to different supra-molecular struc- 
tures, all consisting of arrangements of /5-strands. Apart from single /3-sheets, lam- 
inated multi-sheet structures were found near free-energy minima for the six-chain 
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system. It should be pointed out that the six-chain system is still too small to per- 
mit the formation of, for example, a barrel-type structure. It will therefore be very 
interesting to try to extend these calculations to larger system sizes. 
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